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Abstract. We study discrete-time quantum walks on a half line by means of spectral analysis. Cantero et al. 1 
showed that the CMV matrix, which gives a recurrence relation for the orthogonal Laurent polynomials on the unit circle 
[2] , expresses the dynamics of the quantum walk. Using the CGMV method introduced by them, the name is taken from 
their initials, we obtain the spectral measure for the quantum walk. As a corollary, we give another proof for localization 
of the quantum walk on homogeneous trees shown by Chisaki et al. [3]. 

1 Introduction 

The random walk (RW) is one of the most important models to analyze problems in various fields. The quantum 
walk (QW) is a generalization of the RW [4] with a hope that the QW plays such a role in the quantum field, 
fn fact, it has been shown that there are many results on useful applications of quantum speed up algorithms 
of QWs (see [5] for a nice review) and approximations of QWs describing physical processes given by Dirac and 
Schrodinger equations [fJl [71 15] 15] fTU] . 

The time scaling order of weak convergence of QW with a space-homogeneous quantum coin is the square 
of corresponding classical one. That is, while a random walker spreads in proportion to the square root of 
time, a quantum walker spreads in proportion to time. In general, the shape of time scaled distribution is 
an inverted bell-shape on a bounded support. See for the explicit expressions for the densities, for example, 
[TTjrrj] (one-dimensional infinite lattice Z), [T5] (two-dimensional infinite lattice Z 2 ), and [3] (^-homogeneous 
tree T K ). The review over the detailed discussions on the weak convergence of QWs can be seen in [T3]. Another 
important property of QW is localization which is defined by limsupj^^ P(X t — 0) > in this paper, where X t 
is the QW at time t starting from the origin. The emergence of localization on Z is shown in the case of three- 
and four- dimensional quantum coin [15[I16] by using the spatial Fourier transform, and spatial inhomogeneous 
quantum coin by using the spatial Fourier transform |I7j , and a path counting method |I8j . Through numerical 
simulations, localization occurs with a periodic perturbation [19 . Localization in an aperiodic perturbation 
is shown in [5U] by using a "self duality" property of the QW. Localization on the other graphs have been 
investigated, for example, Z 2 [T3l [21], semi-infinite one-dimensional lattice (Z + = {0, 1,2,... })[22] and T K [3] 
by using the spatial Fourier transform and a path counting. 

In this paper, we analyze the spectral measure and the corresponding Laurent polynomials of the unitary 
operator describing the one-step time evolution of the whole system of the QW. In particular, we focus on 
localization of the QWs. The analysis is based on the CGMV method introduced by Cantero et al Q]. The 
name CGMV is taken from the initial of the authors pQ . Path counting methods [TTJ [T2] or the Fourier analysis 
p~5l [23] are useful for the homogeneous case. However, in general, computations in these methods become 
complicated for the inhomogeneous case, for example, one defect model [18]. On the other hand, the spectral 
theory on the unit circle has a long history and has been well studied (see [21 [25l [26l [27j [28]). The CGMV 
method used here is based on the theory. So it is expected that it would skip some involved procedures in 
path counting methods and the Fourier analysis. Indeed, generating functions of the amplitudes corresponding 
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to the limit measures can be easily obtained by using the Caratheodory functions in the CGMV method, as 
compared with usual two methods [3j[22]. Let p be a probability measure on <9B = {zeC:|;z| = l}, where C 
is the set of complex numbers. L 2 AdJti) denotes the Hilbert space of /x-square-integrable functions on 3D. The 
Laurent polynomials {xi( z )}hLo are orthogonal polynomials on dU> obtained by applying the Gram-Schmidt 
orthonormalization to {1, z, z , z 2 , z~ 2 , . . . } with respect to the following inner product: 

U,9) = \ f{z)g{z)dp(z), 
Jzedo 

where f,g G L^dO). By the recurrence relation between Szego polynomials coming from {l,z, z 2 , ...} on 
9ID [25j and the relation between the Laurent and Szego polynomials, a recurrence relation for the Laurent 
polynomials can be obtained (see [3S1H7] for the detailed discussion). It is shown in [28 that the CMV matrix, 
which is a five-diagonal band representation associated with the Verblunsky parameter (ao, ai, a%, . . . ), gives 
the simplest expression for the recurrence relation of the Laurent polynomials. The CGMV method gives us 
the CMV matrix [2] corresponding the QW. The Fourier transform [23] gives us simpler computational method 
for the QW than the path counting method [TT| [P2] in some cases. However since the Fourier transform maps 
the whole state space of the QW to the quantum coin space, it is hard to recover some aspects of the QWs. On 
the other hand, the spectral analysis for the CMV matrix studied in [TJ [57] systematically gives us a more 
efficient computation of the full system of the QW. In this paper, we show that by using the CGMV method, 
the localization of the QW derives from the point mass of the spectral measure on <9B. Furthermore, we give 
a necessary and sufficient condition for localization and the stationary distribution discussed in [Kl [2TJ [5S] for 
some class of QWs using the CGMV method. As a corollary, we obtain the stationary distribution of QW on 
T K which was shown in [3] using a path counting method and the Fourier transform. 

The rest of this paper is organized as follows. We consider the following two cases: (1) the CMV matrix with 
null-odd Verblunsky parameter, i.e., (ao, ai, a2, . . . ) = (a, 0, a, 0, ... ) and (2) null-even Verblunsky parameter, 
i.e., (ao, at±, a.%, . . . ) = (0, a, 0, a, . . . ). In Section 2, we explain the relations between the CMV matrix and 
the corresponding QWs. That is, we clarify a relation between the CMV matrix with null-odd Verblunsky 
parameter and the QW on Z + with a self loop at the origin (Type I QW), and a relation between the CMV 
matrix with null-even Verblunsky parameter and the QW with a reflection wall studied in [3l 122], In Section 3, 
we present results on localization for Types I and II QWs with a general spatial constant coin. We discuss the 
relation between localizations of both Types of QWs and the eigenvectors at the point mass of spectral measure 
for the corresponding CMV matrix. As a consequence, we give another proof for localization of the QW on T K 
shown in [3J. Finally, in Section 4, we summarize on results. 



2 Quantum walks and CMV matrices 

The CGMV method is a new analytical tool for QWs introduced by Cantero et al. [1] The method is different 
from path counting [3J [TTJ [T2] and the Fourier analysis [21 [HI EH [23] ■ We can get all information on QW from 
the eigen system of the unitary operator W which describes the one-step time-evolution for the system. Cantero 
et al. p] have shown that W can be expressed by a CMV matrix C related to Szego polynomials [2 [28] and 
analyze the spectral measure of the CMV matrix corresponding to W . In this section, we define the two types 
of QWs (called Types I and II QWs) and give the corresponding CMV matrix, respectively. 

Let C( QOl ai,...) be the CMV matrix associated with the Verblunsky parameter {a.Q,ai, . . .). The simplest 
matrix representation for C( a0)Q1 ,...) is given by the following form [21 128]: 



C 



(ao, ai,... ) 



Po 







p ai 
-atoai 

P\U~2 

P1P2 





PoPi 
-aopi 
-aia 2 
-aip 2 










-a 2 a 3 

P3^A 
P3P4 






PlP-A 

-a 2 p3 
-a 3 a 4 
-a 3/ 9 4 








p 4 a 5 
— a 4 a5 








P4P5 

-a 4 P5 



(1) 



where ctj E C satisfies \atj\ < 1 and pj = y/l — \aj\ 2 . At first we present the relation between the QW and the 
CMV matrix. We give definitions of Types I and II QWs, respectively. Let the so called "quantum coin" U be 
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u = 



G U(2), (2) 



denoted by 

crr clr 
crl c L l 

where U(2) is the set of 2 x 2 unitary matrices. We define p = \crr\, &ig(c RR ) = cr R , arg(c^i) = cr^, and 
det(J7) = A, where arg(z) is the argument of z G C. By the unitarity of U, we have p = \cll\, A = e^ aR+aL \ 
and c RL = -AcTr. 

1. Type I QW 

The total space liS 1 ^ for Type I QW is generated by a standard basis {|0, S), 

\0,L), \l,R),\l,L), \2,R), |2,L},. . .}. The time evolution W^' U ^ with quantum coin U is denoted by 

wVW\o t S) = c RR \l,R)+c LR \0,S), W^\0,L) = c RL \l,R) + c LL \0,L), 
wVM\x,R) =c RR \x + l,R) + c LR \x-l,L) (x>l), 
W iLU) \x,L) = c RL \x + l,R)+c LL \x-l,L) (x > 1). 

2. Type II QW 

The total space W 11 ' for Type II QW is generated by a standard basis 

{|0, L), |1, R), |2, R), 1 2, L), . . .}. Let 7 e M, where R is the set of real numbers. The time evolution 

\y(u,u) is denoted by 

W {II ' u) \0,L) = e^\l,R), 

W^ I ^\x,R) = c RR \x + l,R) + c LR \x-l,L) (x>l), 
W^ I ^\x,L) = c RL \x + l,R) + c LL \x-l,L) (x>l). 

In this paper, we restrict initial states of Types I and II QWs to = a|0, S) + /3|0, L) and = e ,;,5 |0, L), 

respectively, where \a\ 2 + |/3| 2 = 1 and S G R. Define a total weight of the passage from position y to position 
x at time t with the quantum coin U by 

^x,'^ (t) = ^ I {(x,d 1 ),(y,d 2 )£-H( J )}(dl,d2) 

di,d 2 &{R,L} 

(x, d x \ (V^)' \y, d 2 )\ IdiXda] (J G {/, //}), 

where = T [l,0], and \L) = T [0, 1] and Ia(x, y) is the indicator function of A, that is, Ia (x, y) = 1 
((x,y) G A), = ((x,y) ^ A). Here T is the transposed operator. Let and Xj ^ be Types I and II 

QWs, respectively. Then the probability that a particle is measured in location x at time t is defined by 
p^ x [J,u) = ^ = pW^Wlp (j e {/,//}), where ip[ J) = T [a,f3] (J = I), = e lS (J = II). To get the matrix 
representation of W^ I ' U \ we give a one-to-one correspondence between the basis of %^ and {0,1,2,...}, 
such that (0,5) o 0, (0,i) O 1, (k,R) o 2fc, (fc,i) 2fc + 1 (fc > 1). Put A, = diag[l, A^ \ A^, . . . } 
with A^? = e~ tka ' R (k > 0), A^_i = e lkc,L (k > 1), where diag[Ao, Ai, . . .] is a diagonal matrix, that is, 
(diag[A , Ai, . . .]) l m = 5i >m Xi. Then we see 

M/^ c/ )=A / T C (ao , , Q2 , ,... ) A / , (3) 

where a,j = aA~^ +1 " 2 (j — even), = (j = odd) with a = c^rA 1 ^ 2 . Equation ([3]) gives a relation between 
the null-odd CMV matrix and Type I QW. On the other hand, for Type II QW, a one-to-one correspondence 
between the basis of U {u) and {0, 1, 2, . . . } is given by (0, L) O 0, (fc, R) O 2fc - 1, (fc, L) O 2fc (fc > 1). Define 
A// = diag[l, Ai /7) , A^ /} , . . . ] with A^ 7) = e -^(^-7) ) xg 1 ^ = ^r-i), Thu s we also see 

^^• c/ ) = e ^A /J C (0 , 6l ,o,6 3 ,...)A IJ) (4) 

where bj = &(e~* 7 A) _ ^' +1 ) (j = odd), = (j = even) with 6 = c/TRAe - * 7 . Equation Q gives a relation 
between the null-even CMV matrix and Type II QW. Denote the Laurent polynomials for C( QOjQ , li ... ) as Uj(z) 
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(j = 0,1,...) and Vj(z) satisfying u(z)Cr a0jCtlt ,,, ) = zu(z) and Ct aoai \v(z) = zv(z) respectively, where 
u(z) = [uo(z),ui(z), . 
is changed to 



{a e lw ,a 1 e 2lw ,a 2 e 3lw , 
as follows: 



and v(z) — T [vo(z), vi(z), . . .]. If the original Verblunsky parameter (ao,cti,a 2 , . 
then the corresponding Laurent polynomials and spectral measure can be rewritten 



u 2k -i(z) -> e ikw u 2k ^(e- ikw z), u 2k (z) -> e -' lkw u 2k {e- lkw z), 
v 2k -i(z) -> e-^uajt-iCe-^z), t* fc (z) -> e^^e^z), 



(5) 
(6) 
(7) 



Let a £ C with lorl < 1 and 



C{a) = 



From the definition of the weight of a passage and Eqs. (1) and (2) for A = 1, 7 = case, we should remark 
that 



B£{?<"»(t) = 



( C («.0,a,0,... j)^ ( C (a,0,a,0,...)) 
C (a,0,a,0,... )) 2l/i2x+1 ( C (a,0,a,0,.. )) 2y+1 , 2a;+ 

( C (0,a,0,a,... )) 2l _ 12rl ( C (0,a,0,a,... )) 2k _ 1>2 



u (0,a,0,a,...) 



L '(0,o ! ,0 : a,... ) 



(s,v>o), 



(z,v>o), 



where if i,j < 0, then we put (c^ ag ai 



...) 



= 0. Thus Eqs. (J3j- (jTj) and orthonormahty of the Laurent 



polynomials give the weight of a passage with a general quantum coin U, ^x,^ [t), as follows: 

= (A 1 / 2 )*^*"^ x I>*(^)S^i, O(o)) -D(0), 



(8) 
(9) 



where V = o R - 7, = (a R - a L )/2, D{0) = diag(e ie / 2 , e~ ie / 2 ), a = cTrA 1 / 2 , and b = cERAe'^ with 
A = det(U). As a summary of this section, Types I and II QWs with quantum coin U defined by Eq. (1) can 
be re-expressed by the CMV matrices with the Verblunsky parameters (a, 0, a, 0, ... ) and (0, b, 0, b, . . . ) under 
the correspondences 

(10) 



ci = cl_rA 1/2 , b = c LR Ae 



respectively. The explicit expressions for S^^"^ and can be obtained by the spectral analysis for 



C(a,0,a,0,...) an d ^(0,6,0, 6,... ) 



>x,y 

in the next section. 



3 Localization and point mass of spectral measure 

In this section, we will consider the relation between the spectral measure of the CMV matrix and localization 
of the corresponding QW. Here we define "localization" as follows: there exists € H^ J ^ such that 



lim 



su J >\(^ J \(w( J ' u A t * { Q J) )\>0, (J e {I, II}). 

•00 \ * 



As we will show, the stationary distributions for Types I and II QWs are described by the Laurent polynomials 
at the point mass of the corresponding spectral measure. Let pS 1 ' and [i^ 11 ^ be the spectral measures for 
C( Q ,o,a,o,... ) an d C(o,&,o,6,...)! respectively. Let {xj(z)}°^ and {Xj( z )}'jLa denote the Laurent polynomials of fi^ 
and n^ 11 ^ satisfying 

C(a,0,a,Q,...)X(z) = zx(z), X(z)C(p ibfi ,b,... ) = Z X{ Z )> i 11 ) 
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where x(z) = T [x Q (z),x 1 (z),x 2 (z),...} and x(z) = [xo(z), Xi(z), Xi{z), ■ ■ ■] with x (z) = xo(z) = 1, respec- 
tively. Therefore 



{p(a,0,a,0,...)j 
(p{0,b,0,b,... )) 



ini 



ini 



1*1=1 



z t xi{z)x m {z)d^ I) (z), 



kl=i 



To get the spectral measure we compute the Caratheodory function F^ J \z) (J R, L}) which are given 

by 



FW{z) = lim ^44, F^ n \z) = lim (1*1 < *)» 



(12) 



where aij(z) and Xj( z ) are the Laurent polynomials whose Verblunsky parameters are — ay when the original 
ones are a-j. Put ={fle [-7r,7r) : linvriF^re 49 ) = 00} for J e Then the spectral measure is 

obtained by 



V J )(e 4e ) = «;( J )(0)^+ J2 rn { J) (e )5(e-e )de, 



2vr 



e ee< J ) 



where 



/ J \e) = \imRc(F iJ \re ie )) , m { Q J \9 ) = lim- — -F^\re l0 °). 

rfl V / rfl 2 



(13) 



(14) 



rtl V ' '/ " " "" rfl 2 

Here Re(z) is the real part of z € C. If w^(9) is i 1 integrable, then by the Riemann-Lebesgue lemma and 
Eqs. ([8} ®, S^Q(t) can be given by, for sufficiently large time step t, 



Bj$(*) ~ (A 1 /2)t e iK^-^)/2 



£ 2fe (e i6) °) 



i(^R-^)/2x 2fe ( e ie o)J 1 ( e i0o) 



e l ^-^)/ 2 x 2fc+ i(e^) £ 2fe+1 (e i( >°)£i(e i( >°) 



e -*(^-7)^ fc _ 1 ( e Wo) 
X2fc(e ieo ) 



. (15) 
(16) 



0o est") 

where — means lim^oo \(A(t)) Lm /(B(t)) Lm \ -> 1 for Z,m € {1,2}. Here A(t) and £(*) are 2 x 2 
matrices. 



4 Main results 

Now we give explicit expressions for the limit measures of Types I and II QWs. 

Theorem 1 Let x[ be Type I QW whose the Verblunsky parameter (a, 0, a, 0, ... ) at time t with the initial 
coin state T [a, 0] starting from the origin. The quantum coin is given by 



U 



with a = cTrA 1 / 2 . Th 



en we nave 



lim P(X^ I} = x) 



t— f 00 

where <fr = (ctr — <tl)/2, and 



Rc(a) 2 
1 - Im(a) 2 



CRR Crl 
CLR c L l 



\ ae iH* + /3e-^ 2 ^(a)\ 2 (l + vj{a)) V f{a) 



-ic/>/2 



^^ { yi-hnW-|Re( a) |}. 



(17) 



Here a R{L ) = arg(c fifi ( LL) ) , A = det(U) and p = y/l - \a\ 2 . 
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From the Eq. (|T7|) . we obtain the necessary and sufficient condition for localization of Type I QW on both 
the quantum coin and the initial state in the following. 

Corollary 1 Localization of Type I QW from the origin with the initial coin state T [a,{3] occurs if and only if 

Re(o) + 0, ae^ 12 + pe-**l 2 Vi{a) ^ 0. 

The first condition depends only on the quantum coin and the second one depends not only on the quantum 
coin but also the initial coin state. 
Next, we give the limit measure for Type II QW. 

Theorem 2 Let .3Q be Type II QW whose the Verblunsky parameter (0,b,0,b,. . . ) at time t with the initial 
state e lS starting from the origin. The quantum coin is given by 

U 

with b = ~CLR/S.e~ 11 . The weight of right moving from the origin is e* 7 (7 6 1). Then we have 
lim P{X [ t H) = x) 

t—*oo 

l + (-lf +t \\M{b)\ 2 :x = 0, 

2 \|A/(6)| 2 (l + l/{^ // )(6)} 2 ){^ / )(6)} 2x ;x>0, [ 1 

where = p/\l + b\ and 

\b\ 2 +Re(b) 



Crr clr 
crl c L l 



M(b) = |l + sgn(|6| 2 + Re(&))} 



(l + bf 



Here A = det(U) and p = y/1 - \b\ 2 . 



From Eq. (ITBl , the necessary and sufficient condition of localization of Type II Q W on the quantum coin 
can be seen as follows. Remark that the following necessary and sufficient condition is independent of the initial 
state e tS . 

Corollary 2 Define D = {(x, y) e M 2 : x 1 + y 2 < 1}. Let the Verbulunsky parameter (0, b,0,b, . . .) of Type II 
QW be b = x + iy with (x, y) G D. Then the necessary and sufficient condition of localization of Type II QW is 
given by 

(x,y)£{(x,y)eD:(x + l/2) 2 +y 2 >(l/2) 2 }. 

The Type II QW for a = 2/k — 1 and 7 = (resp. 7 = tt) corresponds to the discrete-time QW on k- regular 
tree T K with the uniform initial qubit ip$ (Case A) (resp. weighted initial qubit (p$ (Case B)) at the origin [3], 
where tp^ = T [l/^/K,w K /y/n,w 2 /y/K, . . . , w^ -1 / V^l with w K — 1 (J = I), = e 2 W K (j = As a consequence, 
we can give the stationary distribution for the QW on K-regular tree corresponding to the result shown in [3] 
as follows. 

Corollary 3 Let Y t K be the distance from the origin of the QW on T K at time t. Then we have 

1 + (-1)*+- [ C(k) : x = 0, 



lim P{Y t K =x) = — x 

t— >oo 



where 

C(k) = 




Case (A), 
Case (B). 



Proofs of Theorems 1 and 2: We show explicit expressions for the spectral measures and corresponding Laurent 
polynomials for Types I and II QWs in Appendices A and B, respectively. Then substituting Eqs. (TK)1) . (pH]) 
into Eq. (|15p and Eqs. into Eq. (IT51) . implies the desired conclusion with respect to localization of 

Types I and II QWs, respectively. □ 



(> 



5 Summary 



We briefly summarize our results. In this paper we investigated localization and the limit distribution for 
Types I and II QWs by using the CGMV method introduced by Cantero et al. pQ. We explicitly computed 
the spectral measure of the CMV matrix corresponding to the QW. From the point mass of the measure, we 
showed localization of the QW. Furthermore, we obtained the necessary and sufficient condition of localization 
with respect to the quantum coin and the initial state. As a corollary, we gave another proof for localization of 
the QW on homogeneous tree shown in [3] . In addition, we presented the limit distribution for the quantum walk. 
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Appendix A 

We give an explicit expression of spectral measure and corresponding Laurent polynomials for Type I QW. 
The Laurent polynomials {xj(z)}°^ for W^ I,U ^ are given by the following recurrence relation (fTTj) : 



xo(z) = 1, x 2n -i{z) = B+{z)\ n + {z) + B_{z)\ n _{z), x 2n {z) = x 2n -i(l/z), (19) 

where 

±(z- 1 -a)/pT\ T (z) 
X+(z) - X-(z) 

with 

A±(z) = i {z + z- 1 T sgn(Im(e- l l arccos "lz)) v^z - z^ 1 ) 2 + 4|a| 2 } . 

Here Im(z) is the imaginary part of z £ C and sgn(x) = 1 (x > 0), = (x = 0), — 1 (x < 0) for x E K. Then 
we can compute the Caratheodory function by Eqs. (|T2l) and (fT9|) as follows: 

F « W = li m |M = - - s-^-KM a) (2Q) 

j^x 3 {z) ^/{ z - z -i)2 + 4 | a 2| _ 2 Re(a) 
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Therefore from Eqs. (p~3|) and (|20|) . we obtain 



vV-cos 2 ^ 
| sin # + im(a) | 

Furthermore from Eq. ([T4]) . we can calculate Z?^ and m^' as follows: 



= {arcsin(-Im(a))}, mf'ft) = l Re O)l (9 e B^). 



(21) 



These expressions agree with the ones given in [T] . 
Appendix B 

We give an explicit expression of spectral measure and corresponding Laurent polynomials for Type II QW. 
From the reccurence relation (fTTj) . we see that the Laurent polynomials {Xj( z Y\JLo f° r W^ II,U ^ are described as 

Xo(z) = 1, (22) 
X 2n(z) = B^(z)X n + (z)+B^(z)X n _(z), (23) 
X2n+i(z) = B ( °\z)\ n + (z)+B { °\z)\ n _(z), (24) 

where 

± w — \ — — \ — ' ± v z j — r — r\ — \ — n 

A+(z)-A_(z) A + (z)-A_(z) 
Then we compute the Caratheodory function by Eqs. (fT2j) and (|22|) as follows: 



F(-)(,) = (l-6)z-(l-5)*-i (25) 

6z + - v /(z-z- 1 ) 2 +4|6 2 | 



Therefore from Eqs. (|T3j and (|25j) we obtain 



rft s _ V\p\ 2 -cos 2 9 

{0) fcFM J {«e[-^):|c«fl|<p}( Oh 



Furthermore by Eq. (fT4]) , 



- {e e [-7T, tt) : Im ((6 + l)e' e ) - 0} , 
m<, 77) (0o) = M(&)/2 (0 O (26) 



where 



M(6) = |l + sgn(|6| 2 + Re(&))} 



|b| 2 +Re(6) 
(1 + 6) 2 

In particular, in the case of b — ±(—1 + 2/k) which corresponds to the QW on K-regular tree [3J, we obtain 



(27) 



a m h n ^4(«-l)/« 2 -cos 2 r d9 
dn(ff) = (1 - C(k)) 2/wsiiifl I{oe[-,^:\cose\<2V^T/ K }( )^ 



where C(k) = if a = -1 + 2/k, = (« - 2)/(k - 1) if a = 1 - 2/k. 



WW+y-'y (28) 
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